Elliptical Polarization and Probability of Double Ionization 



O 
(N 

1—5 



Oh. 
I . 



o ■ 

Ph. 



in 

o 
o 



Xu Wan£] and J. H. Eberly 
Rochester Theory Center and the Department of Physics & Astronomy 
University of Rochester, Rochester, New York 14627 
(Dated: July 1, 2010) 

The degree of elliptical polarization of intense short laser pulses is shown to be related to the timing 
of strong-field non-sequential double ionization. Higher ellipticity is predicted to force the initiation 
of double ionization into a narrower time window, and this "pins" the ionizing field strength in an 
unexpected way, leading to the first experimentally testable formula for double ionization probability 
as a function of ellipticity. 
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Atoms show an anomalously high degree of electron 
correlation in double ionization when exposed to fem- 
tosecond laser pulses in a wide range of intensities (10^"* 
- 10^^ W/cm^) just below one atomic unit (see recent 
reviews [U, 0). We are concerned here with effects of 
double ionization induced by elliptically polarized light 
(as illustrated in Fig. [1]), over the entire range from lin- 
ear to circular. 

There are two relevant double ionization channels: an 
atom may either lose two electrons one by one, which is 
called sequential double ionization (SDI), or lose the two 
electrons together in an e-e collision between a core elec- 
tron and an already-ionized electron being driven back 
into the core by the reversed phase of the ionizing laser 
field, and this is called nonsequential double ionization 
(NSDI). 

It is considered highly unlikely or impossible for the 
second channel to produce any substantial degree of 
double ionization under elliptically polarized excitation. 
This was confirmed in early experiments Q. It is not 
hard to understand because a collision of the two elec- 
trons will be unlikely or impossible if the returning first- 
ionized electron is steered transversely off course by the 
ellipticity e. 

However, this scenario is at odds with later experimen- 
tal observations. Characteristic NSDI events have been 
observed under circular polarization with the molecules 




FIG. 1: Illustration of an elliptically polarized laser pulse. The 
higher oscillating curve is the electric field along x and the lower 
curve is the electric field along y. The ellipticity is 0.5 here, and 
the full field is: E{t) = Eof{t)[x sin(cjt + ip) + £ycos(ijjt + </>)] for a 
smooth envelope f(t) and random (f>. 



NO and O2 H and with atomic magnesium Ellipti- 
cal polarization has also recently heen predicted to have 
unexpected non-zero effects in SDI .S] , in agreement with 
experiment Q. Evidently, polarization dependence has 
the potential for providing new insights into the complex 
character of two-electron correlation in double ionization 
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We report here "experimental" evidence obtained for 
the dependence of NSDI on the degree of ellipticity of 
the incident pulse obtained via numerical simulations. 
Fig. [2] shows one result, the distribution of initiating 
electric field strengths for successful NSDI events, for 4 
values of e. One sees in the top curve, obtained for linear 
polarization, that the electric field value at which first 
ionization occurs is widely distributed around a broad 
central peak. However, the data changes in a systematic 
way for larger values of the ellipticity. The distribution 
of ionizing fields splits into two peaks as e reaches 0.4 
or 0.5 and becomes quite narrowly localized at less than 
half the distribution's previous peak value in the limit of 
circular polarization. 




FIG. 2: Distribution of field strengths at the time of first ion- 
ization for successful NSDI events, for several different laser field 
ellipticities. 
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We have obtained this data in a series of simulations 
that can be described in the same way as laboratory ex- 
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periments. That is, our simulations represent a "labora- 
tory" in which a laser is focused into a target volume 
with one atom. The peak laser intensity is 6 x 10^^ 
W/cm^. After 10^ laser shots, assuming 100% collec- 
tion efficiency of ions created, there is data representing 
10 million laser-atom exposures. From these, evidence 
for double ionizations must be extracted and analyzed 
in a systematic way. Since polarization dependence is 
of interest, an additional 10 million laser-atom exposures 
must be generated for every important value of ellipticity 
e, say in 10 equal steps between and 1. 

The time at which the first ionization occurs is also 
widely distributed within the laser half-cycle at which the 
peak field reaches its ionizing value, as Fig. |3]shows. And 
the time of first ionization changes in a similar fashion 
to the change in field strength, as e increases, first split- 
ting and then narrowing and ending on a value almost 
exactly half a cycle earlier. This coordination of field 
strength behavior with the timing behavior is not hard 
to understand. The peak in timing should closely corre- 
spond to the peak in the distribution of field strengths, 
and if one is broad or narrow the other should be as well. 
In addition, the jump in timing by half a cycle is eas- 
ily explained, since one expects NSDI events to originate 
near to field peaks, and these occur every half cycle. 
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FIG. 3: Data showing first ionization times for successful NSDI 
events, for several different laser field ellipticities. 

However, the dependence of timing on ellipticity is un- 
expected and we believe not observed or previously sug- 
gested by theoretical considerations. It represents a va- 
riety of NSDI control not previously considered to be 
available. It is known that trajectories with different 
timing can play significantly different roles in high-field 
physical processes such as harmonic generation as well 
as double ionization. So far as we are aware, almost all 
previous investigations of NSDI timings have been re- 
stricted to linear polarization (an exception is the study 
by Shvetsov-Shilovsky, et al. Q). 

To begin to explain this dependence on ellipticity we 
make use of our demonstration that essentially ev- 



ery successful NSDI trajectory has transverse drift in- 
duced by the minor-axis component of polarization that 
is present, but which has been exactly compensated by 
a counter-acting transverse velocity present at initiation. 
Non-compensated trajectories are not able to recoUide 
and are simply absent from any record of NSDI events. 
We believe that the same compensation mechanism is re- 
sponsible for the ellipticity effect shown in Figs. [5] and 

H 

For simplicity, as a first approximation we can assume 
that the transverse momentum distribution available to 
the first electron at its instant of ionization is Gaussian: 
P{vy) exp{—Vy/Avy). The compensation we referred 
to above is accomplished when one of these Vy values 
matches the e-promoted transverse drift velocity, which 
equals sEi/ij, where Ei is the field strength at ionization. 
Thus we will substitute Vy — eEi/uj in the Gaussian ex- 
ponent. 

The overall probability of an NSDI event is then rea- 
sonably estimated as the product of the first electron's 
"release probability" times the Gaussian probability of a 
velocity compensation: 

PNSDiiE,,s) ^ e-Q/^^ X e-(^^^/")VA.2, (i) 

where the dependence of the first factor on ionization 
field strength is chosen in an obvious way, to mimic the 
main effect of Ei in tunneling ionization, while Q is a 
parameter related to ionization potential and varies by 
species. The same is probably true of Aw^. 

To examine in a naive way the effect of e on the Ei 
dependence of ionization we simply plot separately the 
two terms in the exponent in ([T]) as a function of Ei. 
This is done in Fig. 21 where the second exponent con- 
tributes three curves for different values of e. The max- 
imum NSDI probability will come from the minimum 
value of the exponent, and from the graph in Fig. 2] 
it is clear that this minimum occurs at, or very close to, 
the crossing of the curves. That is, a point substantially 
away from a crossing point will find either one or the 
other of the curves rising and making the exponent value 
a lot larger. An additional restriction is that Ei should 
not exceed the critical over-the-barrier field Ec because 
first ionization is expected to saturate under this field 
strength. Ec determines a critical ellipticity value Sc' for 



E =1 




FIG. 4: The two exponents in JTJ are separately plotted as func- 
tions of El . The second exponent is shown three times, for high, 
critical, and low values of ellipticity. The text explains the signifi- 
cance of the curve crossings. 
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e > Ec, the most probable Ei value is determined by the 
crossing point; for e < Ec, the most probable Ei value is 
simply Ec- 

However, the more important feature is the variation in 
the slopes of the curves at the crossing points. The slope 
is much smaller for the lowest values of e, corresponding 
to near-linear polarization. Oppositely, the slopes are 
so high for crossings at or near to values of e « 1 that 
the value of e practically pins the value of Ei . In other 
words, low ellipticities permit a wide range of Ei values 
to create successful NSDI events, whereas near-circular 
ellipticities constrain Ei to a very narrow range. This 
behavior is exactly what is needed to correspond to the 
dependence on e in Figs. [2]and|3l 

The significance of these findings goes beyond the ex- 
planation just discussed. We believe that they expose 
interesting unexplored territory within the NSDI domain 
by showing that new effects appear that depend on el- 
lipticity. Prior to this, in the domain of extensive exper- 
imental activity, mostly confined to near-infrared wave- 
lengths close to 800 nm, few theoretical studies of ellipti- 
cally or circularly polarized pulses have been made (see 

The simulation method that we used to obtain the "ex- 
perimental" data in Figs. [2]and[3]is the classical ensemble 
method that has been described many times (see details 
in [l^). and its validity in interpretation of many dou- 
ble ionization phenomena has been presented [13| and 



critiqued [IJ]. Classical modeling certainly misses true 
quantum features, but since the foundation for visual- 
ization of the NSDI channel is a classical view that at- 
tributes all action immediately after ionization to the 
classical force of the laser field on the freed electron as it 
returns to the vicinity of the ion core [l^, [l^ , a classical 
model is not inappropriate for first analyses. 

In further support, one can say that among theoretical 
approaches the classical modeling used here is the most 
flexible and most widely applicable. For example, it is 
unique up to now among theoretical methods in finding 
agreement with prominently observed features of ion mo- 
mentum spectra in triple ionization [TtI . [18| and with the 
first data on momentum spectra from double ionization 
under elliptical polarization 0j . 

To summarize our use of the method quickly, a micro- 
canonical ensemble of 10^ members is generated using a 
many-pilot-atom method [T^ before turning on the laser 
field. The energy of each 2e member of the ensemble is 
set to be -1.3 a.u., which is close to the binding energies 
of both Xe (-1.23 a.u.) and Kr (-1.41 a.u.) [13, and 
the wavelength is set at 780nm. The width of the fa- 
miliar soft-Coulomb "Rochester potential" 12| is taken 
as a = 1.77, which can be considered the model's single- 
parameter treatment of core effects that are species de- 
pendent (see also (U). Experience with full-dimensioned 
calculations using this method [22j has shown that out-of- 
plane effects can be neglected under current experimental 
conditions and we need only be concerned with the x-y 
plane, as is done here. 
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FIG. 5: NSDI probabilities for 10 values of ellipticity, from our 
numerical experiments. The connecting lines are only to guide the 
eye. An £c value of about 0.5 is estimated from this figure, as 
indicated by the dashed line. 



Our large-ensemble simulations lead directly to NSDI 
probabilities for any value of e, as shown in Fig. [5l The 
results predict a dramatically slower decrease over the full 
range < e < 1 for NSDI rates compared to those indi- 
cated by any previous theoretical considerations known 
to us. As Fig. [S] shows, they fall only 3 orders of mag- 
nitude and remarkable slowly for the higher values of e. 
We believe this provides the first explicitly e-dependent 
explanation for the high NSDI rates in the high-ellipticity 
data from Guo and Gibson and Gillen, et al., mentioned 
already 0,11]. 

Our analysis of ellipticity dependence in NSDI pro- 
duction should be directly testable experimentally. To 
see this we exploit the "pinning" of Ei values mentioned 
above to eliminate Ei in favor of e in the exponent of llj . 
This provides a simple, even simplistic, formula for NSDI 
probability as a function of ellipticity. A quick check 
shows that for values of Ei that are strongly pinned, for 
e > Ec (the value of which can be estimated to be about 
0.5 from Fig. [5]), our expression predicts that the expo- 
nent obeys a power law. Given expression ([1]) as writ- 
ten, the power is -f2/3, and the consequent distribution 
exp{—l3e^^^) is fit in the high-e tail by a /3 value in the 
neighborhood of 4 a.u. 

In conclusion, we have shown that the degree of ellip- 
tical polarization of intense short laser pulses is related 
to the timing of strong-field non-sequential double ion- 
ization. Higher ellipticity is found to force the initiation 
of double ionization into a narrower time window, and 
this in turn "pins" the ionizing field strength in an unex- 
pected way. Among the consequences is an experimen- 
tally testable formula for double ionization probability 
as a function of ellipticity, which predicts a remarkably 
slow decrease in probability at high ellipticities, and we 
believe answers for the first time the standing need for an 
integrated theoretical explanation of the high NSDI rates 



4 



under highly elhptical polarization in well-cited experi- 
ments [3nil- The relationships discovered and reported 
here are generic, and were not tailored to a specific atom. 
But, as we have already mentioned, our classical model- 
ing has been adequate in earlier studies for good semi- 
quantitative correspondence with multi-species NSDI ef- 



fects, as recorded under linear polarization, and we ex- 
pect the same will be true for non-zero ellipticity. 
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